Compact source of narrow-band counterpropagating polarization-entangled photon 
pairs using a single dual-periodically poled crystal 
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We propose a scheme for the generation of counterpropagating polarization-entangled photon pairs 
from a dual-periodically poled crystal. Compared with the usual forward-wave type source, this 
source, in the backward-wave way, has a much narrower bandwidth. With a 2-cm-long bulk crystal, 
the bandwidths of the example sources are estimated to be 3.6 GHz, and the spectral brightnesses are 
more than 100 pairs/ (s GHz mW). Two concurrent quasi-phase-matched spontaneous parametric 
down- conversion processes in a single crystal enable our source to be compact and stable. This 
scheme does not rely on any state projection and applies to both degenerate and non-degenerate 
cases, facilitating applications of the entangled photons. 

PACS numbers: 42.65.Lm, 42.50.Dv,03.67.Bg 



I. INTRODUCTION 

Polarization-entangled photons play a key role not 
only in testing the foundations of quantum mechan- 
ics [I| but also in various photonic quantum technolo- 
gies [^l- A compact, robust, and high-brightness source of 
polarization-entangled photons is therefore desirable for 
practical implementation of a variety of entanglement- 
based applications. 

Spontaneous parametric down-conversion (SPDC) in 
nonlinear crystals is a successful technique to generate 
polarization-entangled photon pairs. A typical method 
involves using the type-II birefringence phase-matching 
(BPM) in a nonlinear crystal jSjj , such as beta barium bo- 
rate (BBO). However, only a small fraction of the total 
emitted photons, the intersecting locations of two non- 
overlapping cones, are polarization-entangled, and there- 
fore, such a source is inefficient. A more efficient source 
consists of two type-I nonlinear crystals via BPM 0, Q , 
from which polarization-entangled photons are emitted 
in a cone. However, generally only a small fraction of the 
cone is collected for use, and thus, such a source is again 
less efficient. 

One way to solve the inefficiency problem in the 
conelike sources is by means of quasi-phase-matching 
(QPM) [101 in periodically poled (PP) crystals [1], such 
as periodically poled lithium niobate (PPLN) and pe- 
riodically poled potassium titanyl phosphate (PPKTP). 
QPM has advantages over BPM due to its higher effi- 
ciency and the fact that it enables flexible frequency- 
tunable processes. In particular, QPM enables the pho- 
ton pairs in a collinear and beam-like configuration. Con- 
sequently, it is possible to make a much bigger fraction of 
the created photons polarization-entangled than conelike 
sources, thus leading to more efficient sources. However 
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a new problem arises, namely the need to spatially sepa- 
rate the collinear photon pairs. 

A simple method to solve this problem is by using 
dichroic mirrors when the photon pairs are generated 
at substantially different frequencies. In this way, sev- 
eral non-degenerate polarization entanglement sources 
have been designed by coherently combi ning two SPDC 
sources at a polarizing beam splitter |9l4l5l|. by manip- 
ulating polarization quouarts IB], by overlapping two 
cascaded PP crystals ll^ . [is'] , or by two cascaded 
or concurrent [H [2l| SPDC processes in a single 
PP crystal. Theses non-degenerate sources have vari- 
ous applications, for instance, in quantum communica- 
tion [2^. However, in many entanglement-based ap- 
plications, for example, in quantum computation [23|, 
frequency-degenerate polarization-entangled photons are 
required. A straightforward way to build degenerate 
entangled sources based on PP crystals is by separat- 
ing collinear orthogonally polarized photon pairs with a 
beam splitter followed by twofold coincidence measure- 
ment as a postselection [24]. However, this method suf- 
fers a 50% loss. A postselection-free method employs 
interferometers to combine two pairs of orthogonally po- 
larized photons [1^, HB] , but such interferometric sources 
(also the non-degenerate sources in Refs. [9-15]) require 
stringent phase control and stabilization. 

Another problem of SPDC sources lies in the broad 
bandwidth determined by the phase-matching condition, 
which is usually on the order of several THz or hundreds 
of GHz. The broadband SPDC source becomes very dim 
in many applications requiring narrow-band photons, 
such as long-distance fiber optical quantum communica- 
tion ('-^GHz ^2?]), strong interaction of the photons with 
atoms and molecules(~MHz 28], and recently relaxed to 
several GHz [2^ [s^l), and interference of independent 
sources without time synchronization (^GHz [Sri). Pas- 
sive filtering is a straightforward way to obtain narrow- 
band sources [13, HH , but it will greatly reduce the gen- 
eration rate. Cavity-enhanced SPDC can provide high- 
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brightness narrow-band photon paris |32h34| . However, 
additional spectral filtering is required to obtain single- 
mode output due to the broad gain bandwidth. 

In this paper, we succeed in solving all the above prob- 
lems by building a compact and narrow-band polariza- 
tion entanglement source based on the backward-wave 
type SPDC in a dual-periodically poled crystal. The 
backward- wave type SPDC |35l - l37l |. has a much narrower 
bandwidth than the forward- wave interaction. The coun- 
terpropagating photon pair generation has also been ex- 
tensively studied in waveguide structures [38l - l4ll | . More- 
over, it not only has the same advantage as the usual 
coUinear, beam-like output SPDC on photon collection 
and overlapping for possible polarization entanglement, 
but it also does not suffer from the problem of spatial sep- 
aration. Our scheme relies on the coherence of two con- 
current backward-wave type SPDC processes in a single 
PP crystal, rather than any interferometer and postse- 
lection. Furthermore, this scheme can work in frequency 
degenerate and non-degenerate cases, for which we de- 
sign two experimentally feasible structures, respectively. 
With a 2-cm-long bulk crystal, the bandwidths of the 
two sources are estimated to be 3.6 GHz, with spectral 
brightnesses of 115 and 154 pairs/(s GHz mW), respec- 
tively. 

The rest of this paper is organized as follows. In 
the next section, we give a description of the dual- 
periodically poled crystal and design the structures re- 
quired in our scheme. In Sec. IIII[ we introduce our 
scheme and make detailed calculations on the sources we 
propose. Sec. lIVI are conclusions. 
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FIG. 1. Schematic of a dual-periodically poled potassium 
titanyl phosphate crystal. Gray and blank areas are inverted 
(— x'^') and background positive (x''^') domains, respectively. 



their Fourier expansions can be written as 
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respectively, where the reciprocals are 

2m7r 2n7r 

— 7 I Ct„ — — , 
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■ sin(m_Di7r), Q„ 



■ sin(ni?27r). (4) 



with Ai and A2 (Ai < A2) denoting the two modulation 
periods, Di and D2 representing the duty cycles, and 
nonzero integers m and n indicating the orders of recip- 
rocals. Then we can write the dual-periodic structure 
as 



g{x) = gi{x)g2{x) = } Qm,ne 



II. DESCRIPTION OF A 
DUAL-PERIODICALLY POLED CRYSTAL 



QPM originates from modulation of the second-order 
nonlinear susceptibility x^"^^ ■ It has been advanced to a 
variety of domain structures which allow multiple and 
flexible nonlinear processes in a single crystal, leading to 
compact and integrated devices. A dual-periodic struc- 
ture is one of the QPM structures, which permits two 
coupled optical parametric interactions [l^, H^. Here, 
taking the potassium titanyl phosphate (KTP) crystal 
for example, we design a dual-periodic structure to sat- 
isfy two concurrent SPDC processes. Hp ^ Hg + Vi, and 
Hp Vs + Hi, where p, s, i represent the pump, signal, 
idler fields, respectively, with H (V) denoting the hori- 
zontal (vertical) polarization. 

The schematic of a dual-periodically poled KTP (DPP- 
KTP) crystal is shown in Fig. [1] in which inverted do- 
mains (with —x^"^^) distribute on a -t-x'^"* background as 
a dual-periodic structure. It is formed by twice-periodic 
modulation of x*'^''- Suppose gi{x) and 32(2;) are two pe- 
riodic functions as the sign of nonlinearity x''^'' ■ Then 



(5) 



where 



^ sin(ml?i7r) sin(rtD27r), (6) 

Gm.n —Gm + G„ = — ! . (7) 

Ai A2 

Then the modulation of the second-order nonlinear sus- 
ceptibility x^^^ can be described as 



X 



(2) 



(x) = dg(x) = dy^g,n,ne 



(8) 



where d is the effective nonlinear coefficient. 

An arbitrary twice-periodic modulation could result in 
smaller domains which may make fabrication more diffi- 
cult. A straightforward way to avoid the unwanted small 
domains, is by designing the structure such that 

'A2/A1 = V2 

Di = l/2 , Z is an integer bigger than 2, (9) 

D2^[l/2\/l 
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where [-J is the floor function to get the integer part of 
a number. In practice, this condition can be satisfied by 
tuning the temperature and wavelengths. 

We consider the pump wave vector along the x direc- 
tion and H {V) in the y (z) directions. By choosing 
the right wavelengths and temperature we are able to 
obtain the following QPM conditions for two backward- 
wave type SPDC processes 

Afci = kp^H - ksM + hy - Gni^^m = 0, (10) 

Afca = kp^H - ks,V + h,H - Gm2,n-i = 0, (11) 

where Gmi,ni and Gm2,n2 are given by Eq. ([7]) in the 
case of {m,n} = {mi,ni} and {m,n} — {m2,n2\, re- 
spectively. Here we require the two SPDC processes to 
have the same signal frequency Us and the same idler fre- 
quency uji, with the energy conservation condition ujp = 
where u)p is the pump frequency. In addition, as 
we shall see in Sec. Illli we require mini = ±m2n2. In 
order to show the experimental feasibility of such a struc- 
ture, in the following we design two possible structures 
based on the temperature-dependent Sellmeier equation 
given by Emanueli and Arie [4J]. 

We first design a structure for a degenerate source of 
\p — 655 nm, Xg = Xi = 1310 nm. Such a source could 
find applications in long-distance fiber-based quantum 
information processing, as the wavelength of the photons 
is in the the second telecom window. At a working tem- 
perature of 75° C, we get the two reciprocals for QPM 
as G3.1 = 17.47 fim^^ and Gs.-i = 18.24 /im^^, corre- 
sponding to the two modulation periods Ai — 1.056 /im 
and A2 = 16.36 /zm, respectively. The ratio of the two 
periods is A2/A1 — 15.5, and thus the duty cycle D2 
should be 15/31. 

We design a second structure for a non-degenerate 
source of Xp = 532 nm, — 807.3 nm, Ai = 1560 
nm. This choice is motivated by the photon source re- 
quirements in real-world quantum networks, for example 
photonic memories in quantum repeaters. The shorter- 
wavelength photon of this source can be used for cou- 
pling and entangling atomic systems, and the other pho- 
ton at 1560 nm can be transmitted over a long dis- 
tance in fiber because its wavelength lies in the low- 
loss transmission window of optical fibers. By choos- 
ing the working temperature as 75.5°C, we obtain the 
two reciprocals for QPM as G3.1 = 14.95 /im~^ and 
Ga.-i = 15.93 /im^^, with the two corresponding mod- 
ulation periods as Ai = 1.220 /im and A2 — 12.82 /im, 
respectively, the ratio of which is A2/A1 — 10.5, and 
therefore the duty cycle D2 = 10/21. 

The above two example structures are both within 
current micron and submicron periodic poling tech- 
niques (45l - l47j . In the following section we shall present 
the SPDC process in the DPPKTP crystal and study the 
performances of the two example sources. 



III. GENERATION OF 
POLARIZATION-ENTANGLED PHOTONS 

We consider a classical pump wave illuminating the 
DPPKTP crystal with a length of L in the x direction 
and the interaction volume denoted by V. The induced 
second-order nonlinear polarization is given by 48] 

Pf ^ (r , t) = eoX?jlE, (r , t)Ek (f, t) , (12) 

(2) 

where en is the vacuum dielectric constant and y;-/, is 
the second-order nonlinear susceptibility tensor, where 
i, j, k refer to the cartesian components of the fields. 
Here, we use the Einstein notation of repeated indices for 
tensor products. The Hamiltonian of the electromagnetic 
system can be expressed as 

H [ d^ff D ■ E + —B ■ b] , (13) 
2 V Mo / 

where fio is the vacuum permeability constant. Since 
D = EqE + P, we obtain the interaction Hamiltonian in 
the parametric down-conversion process 

Hi{t) ^^J^£fP-E 

= £0/ d''rx^^'>Ep{r,t)Es{r,t)E,{r,t), (14) 
Jv 

(2) 

where we replace xljk/'^ with the second-order nonlin- 
ear susceptibility x*-^-* which has the form of Eq. 
for an ideal structure. After quantization of the electro- 
magnetic fields, E{r, t) becomes a Hilbert space operator 
E{r,t), which can be decomposed into its positive and 
negative parts E{f,t) = E^+\r,t) + E^-\r,t). Then we 
can rewrite the interaction Hamiltonianas 

Hi{t) ^so f d'rx^'Hx)E^p+\r,t)Ei-'>{f,t)Et\r,t) 
Jv 

+ H.C., (15) 

where H.c. denotes the Hermitian conjugate part. Here, 
we only write the two terms that lead to energy conserv- 
ing processes, and we neglect the other six terms that 
do not satisfy energy conservation and are therefore of 
no importance in the steady state. Note that neglecting 
these contributions is equivalent to making the rotating- 
wave approximation. 

Since the transverse structure of DPPKTP is homo- 
geneous, we ignore the transverse vectors of interacting 
waves and only consider the interaction along the propa- 
gating direction. We consider the case of signal and idler 
photons in forward and backward directions, respectively. 
Then the negative parts of the field operators of the sig- 
nal and idler Es, Ei are represented by Fourier integrals 
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q=H,V 



(16) 



zj \r ^ 



q=H,V 



where Ej,q — i^Jhwjj (47r£ocrtg(wj)), j = s,i. For sim- 
plicity, here we consider a continuous-wave (cw) plane- 
wave pump with horizontal polarization. In addition, the 
pump field is treated as an undepleted classical wave, and 
thus the positive part of its field operator is replaced with 
its complex amplitude 



jpt) 



(18) 



Then by substituting Eqs. ©,(1161), and (HH]) into 
Eq. we obtain 



hEf 



Aire 
X / dw. 



q=H,V q'=H,V m.n 



nq{uJs)nq, (Wi) 

X g»(^s+'^i-<^p)*g-»(*:a,q-fci,5' -'=p,ff+Gm,„)a; JJ q 

(19) 

For the SPDC process, the interaction is weak, so un- 
der first-order perturbation theory the state evolution 
from time t' to t can be written as 



1 /■* ^ 

l^") = |vac) + — Hi{T)dT\yac) 



(20) 



Considering steady state output we may set t' — — oo 
and t = (X). Then we have 

fire^^s+^.-u^^)r ^ 2tt6{uJs +UJ^~ UJp), (21) 

which gives the energy conservation relation 



(jJs + UJi ~ UJp — 0. 



(22) 



The integral over crystal length can be calculated as 



^2.g-»(fca,,-fc,,,,-fcp,«+G,„,„)x ^ LhiLAkqq,), (23) 



where Akqqi = kp^n — k^.q + ki^qi — Gm.n and the h- 
function has the following form 

1 — c~*^ ■ X X 

h(x) = = e"* 2 sine-. (24) 

^ ' ix 2 ^ ^ 

h(LAkqqi) determines the natural bandwidth of the two- 
photon state, as we shall see. In the case of infinite crystal 
length, Eq. ([25)) becomes a (5-function, thus leading to the 



momentum conservation, i.e., the perfect phase-matching 
condition, Akqq> = 0. 

Suppose that perfect phase matching conditions given 
by Eqs. ([TU| and (ITT]) can be satisfied at frequencies ftg 
and fli, with corresponding wave vectors Ks^h, Ksy^ 
Ki^H, and such that 

n, + n,^up, K.^^^IhlMk^ (25) 



with j = s, i and q = H,V. Due to the existence of 
bandwidth, and constrained by Eq. ([22]), we let 



(26) 



where \v\ <C j ~ s,i. Then in the case of the QPM 
conditions given by Eqs. ([TOl) and ([TT|) . we can write the 
state of SPDC as 

I*) =|vac) + AnvdHvL J duh{LAkHv)a\ ui^s + 

X aj y(rii — v)\va,c) + AvHdvHL j di'h{LAkvH) 
X dl yiVts + v)d\H{^i - 1^)1 vac), (27) 

where 

Ad mi'K 
dnv = dgrm,ni = — :— sin sm(niD27r), (28) 



dvH — dQm2,n2 = 



TT^mini 2 
Ad . m27r 



2, "2 2 



7r^m2n2 2 



sin — — sin(n2£'27r), (29) 



2c V ns.Hn,,v ' 



2c V ns.vriLH ' 



(30) 
(31) 



with Uj^q denoting the refraction index of a photon with 
polarization q at frequency f2j. Here Anvdnv and 
^yffC^vff are slowly varying functions of frequency, which 
have been taken outside the integral. 

We can see that the maximally polarization- 
entangled state can be obtained under the condition 
of \AHvdHvh{LAkHv)\ = \AvHdvHh{LAkvH)\- The 
condition of duv — dyn = d' can be satisfied straight- 
forwardly by choosing mini — ±7742^2. In the following, 
we make calculations on h(LAkHv) and h(LAkvH), i-e., 
the spectrum of the photon pairs. In other words, the 
two-photon correlation time is on the order of several 
hundred picoseconds. 



A. Characterizing the spectrum of photon pairs 
generated from our source 

We first expand the magnitudes of the wave vectors for 
signal and idler photons around the central frequencies 
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ris and fli respectively, up to first order in v, 



l,q 



(32) 
(33) 



where Uq{rtj) — dilj/dKj^q are the group velocities of 
signal and idler photons at central frequencies, with j — 
s, i and q ~ H,V. Therefore we obtain 



AkHv = —i^Shv: Shv 



1 



1 



1 1 



, (34) 
■ (35) 



We thus obtain the joint spectral densities for the com- 
ponents \H,V) and \V,H), 



|MiAfcw)|^-sinc2^^^, 



(36) 
(37) 



and the corresponding bandwidths are Aujhv ~ 
1.77tt/{LShv) and Aojvh ~ l-777r/(LS'yjj), respec- 
tively. Compared with the usual forward-wave type- 
II SPDC under the same conditions on crystal length 
and frequencies [i^, the backward- wave source has a 
much narrower bandwidth, with a reducing factor of 

More explicitly, we consider the two example struc- 
tures given in Sec. [TTl and the crystal length is set to 2 
cm. For the degenerate source, we get the bandwidth 
Aujhv ~ AuJvH ~ 27r X 3.66 GHz and the reducing 
factor is 41. For the non-degenerate source, we obtain 
the two bandwidths as Aujhv ~ 27r x 3.61 GHz and 
AuJvH ~ 2tt X 3.63 GHz, corresponding to reducing fac- 
tors of 25.9 and 78.2, respectively. Note that, compared 
with the asymmetric spectrum in the forward-wave case, 
our backward- wave source has an almost symmetric spec- 
trum. 



B. Quantifying the polarization entanglement 
produced by our source 



its normalization constant is the two-photon generation 
rate, given by 



i?= (*2l*2) =d''L' 



\Ahv\'' I diylhiLAkHvi^^W 



VH\ 



dv\h{LAkvH{v))\ ■ (38) 



Substituting Eqs. and ([57)1 into the above equation, 
we obtain 



R = 2Trd'^L 



\A 



HV\ 



\Avh\^ 



C 



Shv Svh 
Then we can calculate the concurrence 

_\{'i'2Wy®CTy\^2)\ 



(39) 



d'^L^\AHvAvH\ 



R 



dvh*[LAkHv{i'))h{LAkvH{y)) 



+ / duh*{LAkvH{i'))h{LAkHv{y)) 



(40) 



By substituting Eqs. (|24l), ([M]) and ([351) into the above 
equation, we arrive at 



C 



2S'n 



(41) 



where 5min = ^^^{Shv,Svh} and (5„ = 
V ns^HUiy/ [nsyrii^H)- 

For degenerate case, i.e., $7^ = ^i, Shv — Svh, 
(5„ = 1, and thus C = 1, so our source can generate de- 
generate maximal polarization entanglement. This fea- 
ture can also be seen directly from the two-photon term 
of the state given by Eq. (P71) . which shows a maximally 
entangled state in the form of i\HV) + \VH))/V2. 

While for non-degenerate case, i.e., ilg ^ f^i, Shv 7^ 
Svh, <5„ 7^ 1, and therefore C < 1, so the entanglement 
is nonmaximal. However, actually there is not a big dif- 
ference between ^„S'//\/and SvH/^n, so the concurrence 
is very near to 1. Explicitly, let us consider the example 
structure given in Sec. |lTl the concurrence of the entan- 
glement generated from which is found to be as high as 
0.9978. 



To quantify the polarization entanglement produced by 
our source, we employ a commonly used entanglement 
measure, namely, concurrence ^5(1], whose value ranges 
from zero, for a non-entangled state, to one, for a max- 
imally entangled state. For a pure two-qubit state 
expressed in a fixed basis such as {|00), |01), |10), |11)}, 
the concurrence C = \{ip\ay ® ay\4')\, where Uy is the 
second Pauli matrix '^*) in the same basis. For our 
source, we need to treat the two-photon term of the state 
given by Eq. (P7)) . denoted as |^'2)- Note that the state 
\^2) is unnormalized and the reciprocal of the square of 



C. Generation rate of the entangled photon pairs 

The photon pair generation rate can be estimated from 
Eq. dMl), by substituting Eqs. dSJ), and \Ep\'^ = 

2P/ {eoHpcS) into it, where P denotes the pump power 
and S represents the transverse area of the pump beam. 
Therefore we obtain 

^ _ TTd'^LPnsfl, /I 1 \ 

eoUpC^S \ns,Hni.vSHV ns.vrii^nSvH ) 

(42) 
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Let us consider the two specific example sources, for 
which we set P = 1 mW, S — 0.01 mm^, and L — 2 cm. 
The nonhnear coefficient d' is given by Eq. (|28|) . where 
nil = 3, ni = 1, and the effective nonlinear coefficient d, 
stemming from (i24, is found to be 3.9 pm/V. Then we 
find the generation rate of the degenerate source to be 
421 pairs/s, and thus we get the spectral brightness as 
2ttR/Auj w 115 pairs/(s GHz mW). The generation rate 
of the non-degenerate source is found to be 554 pairs/s, 
corresponding to the spectral brightness of 154 pairs/ (s 
GHz mW). We have to emphasize that the experimental 
value of the photon pair rate and the two-photon spec- 
trum will definitely be affected by the poling quality, such 
as the deviations and fluctuation of the poling period and 
duty cycle [5l|. However, the state of the art of the pol- 
ing technique can enable us to engineer a nearly idealized 
poled structure. 

IV. CONCLUSIONS 

In conclusion, we have presented a scheme for build- 
ing polarization-entangled photon pair sources utiliz- 
ing backward-wave type SPDC processes in a dual- 
periodically poled crystal. Our scheme does not rely on 
any state projection and can work in degenerate and non- 
degenerate cases. The backward-wave type SPDC en- 
ables the entangled photon pairs from our source to trans- 
mit in a beam-like way, exhibiting more efficient pho- 
ton collection and mode overlapping. Furthermore, the 
backward-wave type SPDC has a much narrower band- 
width than the usual forward-wave one. In addition, our 
scheme employs two concurrent SPDC processes in a sin- 
gle crystal rather than any interferometer, and therefore 



our source is compact and stable. By proper engineer- 
ing on the domain structure a complete set of Bell states 
can be achieved directly from this DPPKTP crystal [H^l . 
This implies further applications in integrated photonic 
quantum technologies. 

We have designed two possible DPPKTP structures for 
degenerate and non-degenerate sources, respectively. Us- 
ing a 2-cm-long bulk crystal, the bandwidths of the two 
sources were found to be ~ 3.6 GHz with spectral bright- 
nesses of 115 and 154 pairs/(s GHz mW), respectively. 
Our high-spectral-brightness narrow-band sources should 
find applications in large-scale quantum networks and 
other fields requiring narrow-band entangled photons. 
Furthermore, we have also quantified the polarization en- 
tanglement via concurrence and found that the degener- 
ate source can provide maximally polarization-entangled 
photon pairs while the concurrence of the polarization 
entanglement generated from the non-degenerate source 
is as high as 0.9978. Finally, the two structures are both 
within current manufacture technologies, and thus we be- 
lieve our sources can be realized in experiment. We hope 
our approach can stimulate more investigations on appli- 
cations of QPM on photonic quantum technologies. 
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